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ABSTRACT 


In this thesis we have considered the problem of 
designing adaptive power system stabilisers for a single 
machine and a simple multimachine power system • The solution 
proposed is different from those which are available at present 
in the literature. 

The conventional non adaptive PSS work satisfactorily 
so long as the system operating conditions do not change signi- 
ficantly. But in practice, loading conditions and therefore, 
the operating conditions , do change significantly. So, there is 
need for design of adaptive PSS. Some work on adaptive PSS 
design has been reported earlier [28-30]. They all involve 
on-line identification of system parameters. This makes the 
schemes rather complex. The problem of implementation of these 
schemes on actual power systems has not been investigated to a 
satisfactory degree. 

In the method proposed here, a set of PSS is designed. 
Each PSS in the set is designed for a particular operating 
conditions. All the PSS in the set are simultaneously operating, 
in parallel. The control signal for stabilization is obtained 
as a weighted sum of the output of these PSS. Thus, in operation 
only these weights are to be adjusted, in accordance with the 
operating conditions obtained. Proposal have been made for 
determination of the weights and their efficacy established. 


The proposed scheme give performance superior to the conven~ 
tional nonadaptive PSS for changing operating conditions* At 
the same time, it avoids the complexities in the traditional 
adaptive PSS design and is simple to implement. 

The proposed algorithm has been tested in detail, for 
various operating conditions, for single machine infinite bus 
system. Its applicability in multimachine system has been 
confirmed by studying its application to a" typical 3 machine 
9 bus system. 
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CHAPTER 1 


1*1 INTRODUCTION 

During the past two decades there has been 
continuous effort to improve power system transient and 
dynamic stability characteristics. Use of fast acting 
static exciters and high gain AVR has become a universally 
accepted practice. They improve system transient stability. 
However, it has been observed that these excitationcontrol 
systems introduce negative damping, making the system dynami- 
cally unstable, in some cases. The low frequency oscillations 
which are typically in the range of 0.2 to 2,0 HZ are due to 
periodic interchange of mechanical kinetic energy between the 
rotors of power generating stations, caused by differences 
between speeds of machines. If a control signal can be 
applied which reduces the difference in speed, these oscilla- 
tions can be damped out. So a feedback of speed error can 
act to control the speed and provide damping. From the above 
argument, speed has become the most common and popular stabi- 
lizing signal. Acceleration, accelerating power and frequency 
are the other signals which can be used for the same purpose. 
Comprehensive study of use of all the popular stabilizing 
signals are available in the literature [l,47j* 

In the foregoing, it has been noted that a stabili- 
zing signal is required for the excitation control loop. The 
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auxilliary controller, known by the name power system stabiliser 
(PSS), receives the stabilizing signal from rotor speed, accele- 
rating power or frequency and provides corrective signals at 
the input of the excitation system in order to damp out the 
oscillations in the system. Since we are concerned with 
system response to small scale disturbance, PSS can be designed 
using the linearized model of the power system about an opera- 
ting point. 

1.2 PSS DESIGN OBJECTIVE 

As already mentioned, the reason for instability in 
power systems is the negative damping introduced by high gain 
excitation control system for improvement of transient stability'" 
limit and system voltage regulation. Since we deal with a 
linearised system, where negative damping can be viewed as 
being due to positive real part of a complex eigenvalue, the 
objective of PSS design will be to shift the unstable complex 
eigenvalues to appropriate locations in the left half plane. 

1.3 SYSTEM MODEL 

The model considered consist of a single machine 
connected to an infinite bus through a transmission line, 

( Fig. 2.1). Development of state equations for this model 
is given in chapter 2 and Appendix A-1. We also consider in 
this thesis a 3 machine, 9 bus system shown in Fig ,A 4.1.State 
equations for this system are developed in chapter 4, 
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1.4 OPERATING CONDITIONS 


The real and reactive power outputs and the terminal 
voltages of the generators are taken as operating conditions 
for the design of PSS. As the load changes the real and 
reactive power generation changes, so the PSS designed also 
should change. In the ideal case for every operating point 
there must be a PSS to provide proper stabilization. 

1.5 DESIGN METHODOLOGIES 
1.5.1 Non Adaptive 

In this design methodology the PSS structure and 
parameters, once designed, do not change with Change in 
operating conditions or in system configuration. Such a PSS 
gives good results only if the system operating condition do 
not change much from that used for PSS design. A rather 
extensive list of methods falling into this category is 
available in the literature. Here only a few of them are 
discussed. 


The various methodsused for PSS design can be 
broadly classified into two categories, l) those based on 
classical control theory and 2) those based on modern control 
theory. 

The problem of low frequency undamped oscillations 
was investigated by De Mello and Concordia [2] through the 
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concept of synchronizing and damping torques produced in the 
system at the oscillation frequencies. One is proportional 
to the difference in rotor angle of a machine to that of the 
system and is called the synchronizing torque; the other is 
proportional to the speed of a machine around system speed and 
is called the damping torque. In Ref. [2] the authors used 
the single machine infinite bus model to explain the effect 
of the excitation system on the rotor oscillations, A design 
method for PSS using frequency response technique was suggested, 
to improve damping torque without affecting synchronizing torque 
Bollinger and Laha [s] used the root locus method of design to 
shift the rotor oscillation eigenvalues to stable locations. 
Multivariable frequency response plots were used by Hughes 
and Hamdan [4,5j for analysis and design of PSS, These methods ; 
help to determine the nature of phase lead compensation needed 
for stabilisation purposes. However these methods are appli- 
cable only to single machine equivalent models of the system and 
are not general enough to be used for the multimachine case. 

Among the methods based on modern control theory 
there are mainly two approaches available, 1) those based on 
linear optimal regulator theory and 2) those based on pole - 
placement techniques. 

In the approaches based on linear regulator theory 
[6-8], the stabilizer is designed as an optimal state feedback ; 
controller which minimizes a quadratic performance index. : 
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The resulting controllers in these methods [6-8] improve 
the closed loop system performance, but the inaccessibility 
of all the states for feedback creates a problem in practical 
applications. Suboptimal controllers with the system repre- 
sented by lower order models with only measurable states as 
state variables have been developed [9j. The resulting contro- 
llers are optimal only with respect to the model considered. 
Davison and Rao [lOj developed a design method for PSS using 
output feedback by minimizing an average performance index as 
proposed by Athans and Levine [ll]. De Sarkar and D, Rao [l2] 
designed excitation control for stabilizing a synchronous 
machine through output feedback by first considering state 
feedback obtained through application of optimal control 
theory and then getting a ’minimum norm’ suboptimal control 
law for output feedback . The weighing matrices and the 
output variables are chosen to improve both small and large 
signal machine responses. 

Lately, there has been interest in the application 
of pole assignment techniques to the problem of power system 
stabilization [l3-20]. This approach provides a design method 
for the PSS with the control objective of assigning the closed 
loop eigenvalues to specified stable locations; Among this 
class of methods, there are two basic divisions namely, 
i) methods using state feedback and ii) methods using output 
f eedf ack. 
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The basic principle of pole assignment techniques 
for PSS is shifting the open loop eigenvalues to desired 
locations by auxiliary feedback* If the system is controllable, 
all the open loop eigenvalues can be shifted to arbitrary loca- 
tions with state feedback provided complex eigenvalues exist in 
conjugate pairs. If the system is not fully controllable, the 
eigenvalues associated with the controllable modes can be 
shifted arbitrarily if the complex eigenvalues exist in conj- 
ugate pairs, by state feedback* In the literature there are 
many pole. assignment techniques available with state feedback 
[l3-16j* Out of these. Ref, [l4,15j offer simple technique 
using modal control theory. However, unavailability of all 
the states for feedback always poses a practical problem for 
these methods* 

From a practical viewpoint pole assignment using 
output feedback is more important. If the system is completely 
controllable and observable then all its modes can be shifted 
by output feedback provided the complex eigenvalues exist in 
conjugate pair [21], Most of the literature available has 
considered PSS design as a problem of designing dynamic 
output feedback [22-27] compensator. The number of poles 
that can be assigned depends upon the order of PSS. Usually 
a second order PSS is designed to assign the dominant 
eigenvalues of the system. 

From the literature [l7-20] it is observed that 
min, [n,m+l-l ] poles of a completely controllable and observable 
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system can be assigned to desired location by gain output 
feedback , where .n represents the system order, m and 1 
represents number of inputs and outputs of the system respec- 
tively* Padiyar et al. [22] have presented three methods for 
the design of PSS by pole assignment with output feedback. The 
first method is based on modal control theory. The second 
and third methods involve static optimization for exact assign- 
ment of poles and assignment into a sector, Arcidiacono et al, 
[23] proposed a technique for PSS design which utilizes the 
residues of open-loop transfer function to get a specified set 
of closed-loop .eigenvalues, Fleming et al, [24] designed a 
sequential algorithm for tuning the parameters of PSS with fixed 
poles in' multiraachine power system. Sequential addition of 
stabilizers disturbs previously assigned eigenvalues and hence 
may not give satisfactory performance in many cases. Lefebvre 
[25] proposed an algorithm for simultaneously tuning of all 
the PSS in a multimachine system, Abe and Doi [26,27] developed 
a design technique by combining the frequency response and pole 
assignment methods. Munro and Hirbod [31] presented a very 
systematic design of full rank output feedback controllers 
which allow arbitrary pole assignment in the closed-loop system, 
Lefebvre [33] has recently used the tuning of PSS in three 
machine system. Here, firstly the design are done for comple- 
tely decoupled system, then the PSS parameters are updated as 
the interconnection is slowly introduced. 
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1.5,2 Adaptive PSS Design 

So far, the PSS design techniques discussed were 
for design at a particular operating point to provide the 
necessary damping torque. However, the operating point is 
not fixed. It keeps changing with change in loading condi- 
tions and the performance of the stabilizers, under changed 
operating condition are not always satisfactory. This is 
why the need for adaptation has been felt. Adaptive stabi- 
lizer design to enhance both dynamic and transient stability 
has been developed by Ghosh. [28], Here the controllers deve- 
loped are based on recursive parameter estimation. For this 
the recursive least squares identification scheme was used. 

Also the proposed adaptive stabilizers have been implemented 
on a 16 bit microprocessor. The availability of low cost 
microprocessor chips has made it possible to design adaptive 
PSS. The use of digital adaptive regulators, called self tuning 
regulators, for power system stabilization has been reported 
in the literature [29,30], 

1.6 SCOPE FOR THE PRESENT WORK 

It has been felt that non-adaptive PSS designs are 
not satisfactory under changing operating conditions. Also a 
fully adaptive design involves complexities in implementation. 
The design method proposed in the present thesis gives good 
performance under changing operating conditions. At the same 
time, it avoids the complexities involved in the traditional 
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adaptive design.. The scheme proposed is general enough to 
be applicable for changing system configuration, (e.g., strong 
and weak systems). So , the designed controller will give satis- 
factory performance over a wide range of operating conditions 
and network parameters. 

A chapterwise brief summary of the work done is 
presented below; 

Chapter 2 considers a conventional non-adaptive 
PSS design using dynamic output feedback. The algorithm 
used was developed by Munro and Hirbod [31 J. A single mach- 
ine infinite bus system model is developed. The stabilization 
singal used is derived from angular velocity of generator 
rotor. Since we prefer to go for complete pole placement, the 
minimum order of PSS required is three. Design of such a PSS 
has been outlined, and its transfer function obtained. 

Chapter 3 contains the proposed adaptive PSS design. 
Numerical results are obtained for a single machine infinite 
bus model, at different operating points. The time responses 
show the superiority of the proposed scheme over the conventi- 
onal nonadaptive designs. Quality of performance of the 
proposed PSS with change in system strength is discussed. 

In chapter 4 we have dealt with the adaptive PSS 
design for a multimachine system. A 3 machine, 9 bus system 
model is considered for the purpose. The proposed adaptive 
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PSS scheme is generalised enough from the earlier formula- 
tion for application in the multimachine case. The design 
of PSS for each machine, for a given operating condition in 
the system is done by following the procedure given in Ref, [ij. 
Numerical results are given at a few operating points, which 
indicates the advantage of the proposed scheme over the exis- 
ting ones. 

Chapter 5 containes the general conclusion. 
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CHAPTER 2 

POWER SYSTEM STABILIZER DESIGN 

ft 

2.1 INTRODUCTION 

The important role played by PSS on overall system 
performance has already been stated in the previous chapter. 

In fact since 1966 a number of studies [35-39] clearly indica- 
tegT the beneficial effects of PSS on damping out oscillations 
of small magnitude in systems provided with modern high gain 
and fast acting voltage regulators. In this chaper we demon- 
strate a simple PSS design using a linearized state-space 
model of a synchronous generator connected to an infinite 
bus through a transmission line. The power system model used 
here is taken from Ref. [34]. An algorithm is presented for 
design of PSS based on pole placement technique. Numerical 
results are also presented for a PSS designed at a particular 
operating point. 

2.2 POWER SYSTEM MODEL 

In the most simple case, the power system consists 
of a synchronous generator connected to an infinite bus through 
a transmission line. The system alongwith exciter-voltage 
regulator is shown in Fig; 2.1. The nonlinear system dynamics 
and linearized system equations derived from the nonlinear 
equations about an operating point are given in Appendix A-1. 
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F2&. 2-1 System diagra m 



2.i(t>>S\r\q[e phase equivalent circuit 
of the model 





fi&.t-tcc) Block diagram of excitation 
system 
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The linearized system equations for a typical operating 
point are given below. Derivation of this model and the 
system, data used for this purpose are given in Appendix A-2. 


The state and output equations are, respectively 


Where 


X 


Ax 4-Bv 


y = Cx 

V = [ A u ] 


( 2 . 1 ) 


( 2 . 2 ) 


A= 


and 


B= 


-0,16312 

-0.3278 

0.0 

-0.3' 

0.0 


0.0 

1.0 

0.0 

-26.34 


-47.5885 

0.0 

0.0 

719.0961 

-773.2739 

0.0 

-100.0 

'"0.34722' 





0.0 





0.0 




* • 

0.0 






. (2.3) 


(2.4) 


For the present purpose of this work the feedback variable 
considered is angular velocity (Aco). Thus . 

y= [Am] 

and = [0.0,0. 0,1. 0,0.0] 


....♦(2.5) 


For the above mentioned case, complete pole placement is done. 
Eigenvalues of the system without PSS are given in Table (2.1) 
From the Table it is clear that the system is unstable. Our 
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prime objective in designing PSS is to stabilize the system 
by shifting the unstable dominant eigenvalues to a suitable 
location in the left half plane. 


SI . No . 

Eigenvalues 

1. 

0.533+j 7.107182 

2. 

0.533-j 7.107182 

3. 

-3,734817 

4. 

-97.4773 


Table 2.1 Eigenvalues of the power system 
without PSS. 

2.3 PROBLEM FOPIMULATION 

Let q be the number of poles to be placed and m 
and 1 be the number of inputs and outputs respectively. For 
our model given by eqns. (2,1“2.5) ,n> m+1-1, where n is the 
system order. The problem is to find an output feedback Iqw 

U(s) = -F(s).Y(s) ....(2.6) 

Where F(s) is ,a proper rational polynomial matrix of minimal 
order r, such that all the (n+r) poles of the compo site system 
get placed at specified locations. 

2.4 ALGORITHM AND RESULTS 

Since in this study we use only one output as 
feedback variable, the transfer function- matrix F(s) of the 
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compensator in equation (2.6) is a scalar function. For 
compensator order of r, its transfer function can be repre- 
sented as : 


F(s) 


r r-^1 r— 2 

S 4‘CC‘j S +CCr)S -4“ • • « • « 

12 r 


... (2.7) 


Since we opt for complete pole placement, order of the 
compensator required is r=3 ( Justification of this statement 
can be found in Ref. [45]. So in our case 


F(s) = 


■"3 

s +a 2 S-HX 2 


... ( 2 . 8 ) 


The algorithm used for determining F(s) was develo- 
ped by Munro and Hirbod [31 ] and is presented in Appendix A-3» 
No restrictions are placed on the parameters of F(s) in their 
determination. Table 2.2 gives the desired and assigned 
locations of closed loop poles as well as the parameters of 
the dynamic compensator which accomplishes the assignment. 


2,5 DISCUSSION 

It is observed from the results that all the pole 
locations in the closed loop system are in the left half of 
the s-plane. It is further observed that due to slight 
numerical inacuracies in computation, the pole locations are 
not ei^act, though there is no substantial deviation between 
the desired locations and the actual locations where the 
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Feedback variable 

w only ■ 

No. of poles to 
be placed 

7 Note 

n=4 

r=3 

n+r=7 

Desired location of 

-1.0+j 7.107182 

closed loop poles 

-l.O-^- 7.107182 


-97.47773 


-3.734817 


-20.0, -15.0, -25.0 

Order of compensator 

3 

Compensator T.F, 

-3 . 7135s^-392..5681s^-3155. 578S-6527 . 1 
s^+62.31s^+1290.97s+7573.24 

Locations of assigned 

- 1.0+j 7.080849 

closed loop poles. 

- 1.0- j 7.080849 
-97.47779 

- 3.731152 

-15.1418+j 5.560658 
-15.1418-j 5.560658 
-28.87101 


Table 2.2 


Numerical results of complete pole placement 








Oynamic compensator 
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poles got placed. The traditional design technique uses a lead/ 
lag-lead type of compensator [46j and conventionally the 
order is kept at 2 (Fig. 2.2). Stabilizing signals used are 
accelerating power, frequency and velocity. Here we choose 
only velocity feedback and we observed that a 3rd order 
compensator is better than the conventional 2nd order 
compensator, .so we opted for 3rd order PSS. The addition of 
a third order compensator has increased the number of poles 
of the closed loop system by three. It is difficult to deter- 
mine what should be the values specified for these three extra 
poles. They are specified to be real and negative with large 
magnitudes. The exact locations are specified rather arbitra- 
rily. Two of these three poles have assumed complex conjugate 
locations and the remaining is real. All are stable. The 
algorithm is computationally efficient, also it determines the 
minimal order compensator, which is an advantage. In the 
literature, there are several other methods available [40-42] 
for the same purpose. However, we have chosen the simplest one. 
The same algorithm will be used in the next chapter for the 
purpose of developing an adaptive PSS for a wide range of 
operating conditions,. 
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CHAPTER 3 

ADAPTIVE PSS DESIGN FOR WIDE RANGE OF OPERATING CONDITIONS 
3.1 INTRODUCTION 

The PSS described in the previous chapter 
works satisfactorily so long as the system operating conditions 
do not change significantly. However, the system operating 
conditions do change significantly. It is shown in [2j that 
excitation system has conflicting effects on the synchronizing 
and damping torques as the operating conditions change. Under 
heavy loading conditions , the effect of excitation system is 
to increase the synchronizing torque in the system and reduce 
the damping torque. The reduction in damping torques, parti- 
cularly in systems with long transmission lines, results in 
extremely oscillatory behaviour of the system • There is, 
thus, need for developing compensators which operate reaso- 
nably well over a range of operating conditions. Indeed 
the compensator should function well even when the system 
configuration changes due to faults or any other reason. One 
solution to this problem is a fully adaptable compensator, 
which, on identification of any change in the operating condi- 
tions as well as system configuration, changes the structure 
as well as parameters of the PSS, such adaptation processes 
would require computers for implementation. They should be 
reasonably fast for on line implementation. This method 
however, suffers from the drawback of high implementation 
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cost and complexity of system. So the search for simple 
adaptive schemes continues. In this work we describe one 
such method which might be a possible answer to the problem. 

It is to be kept in mind that more work needs to be done 
before it can be presented as a practical solution to the 
industry. 

3.2 PROPOSITION 

The method we present here is surprisingly simple 
from its physical point of view. Instead of having a single 
PSS, we have several PSS, each designed for a particular 
operating point. All the PSS are active simultaneously and 
each contributes to the AVR-exciter block through a weighting 
factor. Of this whole scheme, only the weighting factors are 
subject to change as the operating conditions change, A block 
diagram representation of the scheme is given in Fig. 3.1, 
which uses five individual PSS simultaneously, to generate 
a composite PSS, 

3.3 IMPLEMENTATION 

To carry out the work proposed in the previous 
section let us consider the space of operating conditions, 
with its two axes as real and reactive powers generated by 
the synchronous generator at any time. Since real and reactive 
powers supplied are governed by the load condition, the opera- 
ting point moves in this space. We select five typical 
operating points placed symmetrically over the full load 
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Fig. 3.1 An Adaptive Composite PSS Scheme 










Reactive Power (P.U.) 
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Fig, 3o2 Operating Space With Five 
Design Points. 
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curve which covers a range of power factors from 0.6 (lag)to 
0,6 (lead), Fig, (3.2), It is unlikely that a practical 
system will be required to operate outside this range. Now 
for each of the five operating points considered, we assume 
that there is only one PSS. We design the PSS using the 
algorithm presented in the last chapter, so that all the 
closed loop poles are assigned fairly close to the desired 
locations in the left half plane. Thus we design five PSS 
for the five distinct operating points considered. Note that 
the system matrix [a] changes with change in real power (P) 
and reactive power (Q) flow . However, since the number of 
inputs to the PSS ( in this case w only) and the number of 
poles to be assigned does not change, the structure of the 
compensator remains unchanged and only its parameters change. 
The transfer functions for these five PSS are as shown in 
Table 3.1. At any instant the weighted sum of the outputs of 
these five PSS is fed into the AVR-exciter block. A simple 
procedure to determine the weightages is presented in the next 
section. 

3.4 WEIGHTAGES FOR PSS 

It is now necessary to develop a scheme for calcu- 
lating the weightages to be used in the composite PSS shown 
in Fig. 3.1. The scheme should be fast and effective, A 
simple scheme for this purpose is proposed below. 
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Table 3.1 Transfer functions of five PSS. 
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In the design space i.e. the space of operating 
conditions, there are five points at which we have already 
designed PSS. Let us call those five points as design points. 
Now while in operation, we locate the operating point at 
any time in the design space and calculate its linear distances 
from the five design points. Then we assign the individual 
PSS weightages in the ratio of the reciprocals of the distances 
of the operating points for which these PSS are designed from 
the operating point currently obtained , such that the sum of 
the five weightages is unity. Thus if the new operating point 
is such that it is at distances d|^,d 2 , ...., d^ from the 
operating points corresponding to the five stabilizers, 

PSSl, PSS2, ..., PSS5 respectively, then the respective 
weightages ...,w^ are such that 


Wi:w2:w3:w4:w5 = j- 



5 

Furthermore, w. >0, i=l,2, ..o,,5 and £ w. = 1 

^ ' i=l ^ 


(3.1) 


The logic behind such a procedure is simple. Thus, further 
away a new operating point is from a particular design point 
in the design space, less is the contribution of that design 


point towards determining parameters of the composite PSS 


for the operating point. 


Numerical Exami 


The five design points can be written, in terms 
of coordinates (P,Q), as 
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Sj^=(l»0,0.0) j S2=(0.8,0.6), S3=(0.,6,0.8) , S4=(0*8,-0.6) , 
S5=(0.6,i-0.8), ....( 3.2) 

Let the new operating point be located at (0.4, 0.7). The 
linear distances of the five design points from this new point 
in the design space are 


dj^ =0.9219544, 63=0, 4123106, d3=0. 2236068, 64=!, 360147, 

d^=l. 513275. ....( 3.3) 

So, weightages given to the design points are in the ratio 


Wi;w2JW3 

= 1.085; 


:w4:w^; 


2.425; 


= .- 1 _ 

*^1 ^2 *^3 *^4 

4.472;0.735;0,661. 



(3.4) 


and W2+W24-W3+W4 +w^ = 1 u .. . 

Thus the weightages for the PSS can be determined to be 
Wj_=0,1157, W2=0.2586, W3=0.4768, W4=0.0784, W5=0.0705 .... (3.5) 
3.5 COMMENTS 

The algorithm used for calculating weightages lacks 
sophistication in the sense that all it cares about is the 
linear distance in the P-Q plane. It does not discriminate 
between a deviation in real power and a deviation in reactive 
power. However, the simplicity of the procedure and negligibly 
small computing time required for calculating the weightages 
is a notable benefit. We have, of course, to check that 
this scheme gives satisfactory results. The fact that we have 
taken the sum of five weightages to be equal to unity is 
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justified because if by any chance, the new operating point 
happens to coincide with any one of the design points then we 
prefer to use only the PSS designed for that design point. Note 
that in such a case the weight corresponding to that design point 
will be unity, while all other weightages become zero. Also 
we have assumed that the weightages cannot be negative. 

3.6 A PROBLEM ENCOUNTERED 

Apparently the proposition made for having five 
PSS is a very simple one and it should not create any difficulty 
in designing. However, while working with this idea,, we 
encountered what may be called a dimensionality problem. This 
is explained below. 

We have considered a simplified linearized model of 
the power system. This model is of fourth order, and includes 
the AVR-exciter assembly. Also we have five designed PSS, each 
of which is of third order and each having distinctly different 
ploes and zeros. So, in general, at any operating condition 
the order of the combined closed loop system will be 4+3x5 =19, 
which is rather high. Analysis of the combined closed loop 
system is necessary to ensure that the system performance is 
satisfactory, and we would like to avoid problems of high order 
system if we can. So at this stage we decided to reduce the 
order of the closed loop system by reducing the order of the 
composilP PSS even though this implied sacrificing some of 
the design freedom . These considerations become even more 
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important in multimachine system PSS design using the present 
approach. The technique used to tackle this problem is illus- 
trated in the next section, 

3.7 A SOLUTION TO THE PROBLEM OF DIMENSIONALITY 

In order to reduce the order of the closed loop 
system we decided to specify identical denominators for the 
five PSS, so that all the five compensator poles are at the 
same locations. This would imply that we /have to do partial 
pole placement in the design. Though in our previous design 
technique both the parameters of numerator and denominator 
polynomials were determined for acheiving the assignment of the 
closed loop poles, it has been observed [43-44] that the contri 
button of compensator poles to the final assignment of the 
closed loop poles of the system is not as significant as that 
due to compensator zeros. Thus the design freedom sacrificed 
is not much. So, we adopted the technique and prespecified 
all the denominator parameters of the compensator. For doing 
so a minor change was necessary in the program which was 
used for earlier design. Following this new algorithm 

we designed five PSS for the five design points and repeated 
the procedure described earlier in this chapter to decide about 
the weightages. The new PSS transfer functions for the five 
design points considered are given in Table 3,2* Here, as it 
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Table 3.2 Transfer functions of five PSS with 
numerical values of parameters. 
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is obvious from Table 3.2, we have decided to retain the poles 
of PSSland fixed the pole locations of all other PSS to 
coincide with those of PSS 1. PSS 1 is designed using complete 
design freedom. This way of specifying the denominators of 
the PSS transfer functions is convenient and has a rationale. 

In this way our aim to reduce the order of the closed 
loop system has been acheived, Now, for any new operating point, 
we have the closed loop system of order 4+3 = 7 which is signi- 
ficantly lower than the order of 19 in the case of full design 
freedom. 

3*8 NUMERICAL RESULTS 

We have tested the proposed algorithm at 25 test points 
all over the design space, Fig. 3.3 . The results obtained at 
every point were satisfactory. In Table 3.3 we show some of 
the typical results obtained while testing the program, based 
on this algorithm, on a DEC- 1090 computer system. 

From the results in the tabular form, the advantage 
of having the composite PSS comprising of the five individual 
PSS is not very apperant. To realize how superior the new 
algorithm is, we have to analyse the closed loop system with 
only one PSS ( say PSS l) at any test point in comparison with 
the composite PSS at the same test point. We performed this 
analysis with 3 kinds of arrangements: 
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i) Only PSS 1 operating, 

ii) Only PSS 2 operating. 

iii) Composite PSS operating. 

The results are given for four different operating 
conditions, in the form of time response of the closed loop 
system, when a disturbance of fixed magnitude is applied at 
time instant of zero* The state variables d (angular position 
and (A) (angular velosity) are chosen for displaying the system 
response [Fig. 3,4-3.11]. 

3.9 STRONG AND WEAK SYSTEMS WITH THE COMPOSITE PSS 

To study the system with change in system 
configuration, two cases are considered, namely, strong 
system and weak system. The results produced above were for 
a weak system, with transmission line reactance of 0.7 P.U. 

The system designed has been analysed for strong system also, 
with Xj = 0.2 P.U, Improved system performance was observed. 
It has also been observed that even in the case of a very weak 
link, with = 1.4 P.U. , the system remained stable and 
the performance was satisfactory. Table 3,4 justifies these 
statements. 

3.10 COMMENTS AND DISCUSSION 

From the time responses drawn, the superiority of 
the proposed scheme of composite PSS over the conventional 
scheme of a single PSS is clearly evident. With change in 
system operating conditions, disturbances in the system due 
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Fig. 3.4 Step response for rotor angle (A6). Real power = 0.6 PU, Reactive power = -0.8PU 
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Fig. 3.! 
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Fig, 3.6 Step response for rotor angle (A6). Real power =0.4 PU, Reactive power = 0.7 PU 
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Fig. 3.7 
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Fig. 3.9 Step response for angular velocity (Aw), Real power* 0.6 PU, Reactive power* 0.8 PU. 
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Rg. 3.10 St«p response for rotor angle (A6). Real power *0.4 PU, Reactive power* -0.3 Pi 
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Fig. 3.1! Step response for angular velocity (Ato), Real power = 0.4 PU, Reactive power* -0.3 PU 
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to changeover to a different set of weights will not be much 
It is ofcourse, possible to improve on the scheme proposed 
for selecting weightages or develop alternate schemes. 

By prespecifying the poles of the compensator we 
have lost some of the design freedom and thereby the final 
locations of closed loop poles are affected slightly. Howevei 
as it has been observed in Ref. [43], this kind of compensa- 
tion is not very sensitive to the compensator pole locations. 

On the other hand, by doing so, we have acheived a great ad- 
vantage of always having a stable compensator. If the denomi- 
nator parameters of the compensator are free to assume any 
value then in some case it might generate a compensator with 
unstable poles. A stable system made by an unstable corapens- ; 
ator [ 44 J is undesirable. The question of the best choice of ; 
the composite compensator poles is a relevant and interesting : 
problem. This, however, has not been investigated here. 

Now questions may arise as to why we have chosen 

I 

five PSS to make the composite PSS? What happens if more 
number of PSS are employed? Is there any limit on the number i 
of PSS? The answer to these questions are not easy at the | 
present moment; however »it may be attempted in the following 
way : i 

The performance of the system will no doubt improve | 
with increase in the number of compensators. This is because 
with increase in the number of design points, the design space 
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is covered in greater detail. Also the choice of more 
appropriate design points will be an added factor which 
can be investigated. But the question of upper limit on the 
number of PSS must be approached from a more practical point 
of view. It makes sense to have as few PSS as possible in the 
composite PSS from the point of view of PSS hardware and/or 
software, as the case may be. We have found that five 
individual PSS give satisfactory system performance over a 
wide range of operating conditions. It is felt that the 
controller complexity has not been inordinately increased 
by this number of PSS, 
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CHAPTER 4 

ADAPTIVE PSS FOR MULTIMACHINE SYSTEM 

4 • 1 introduction 


So far we have applied our proposed algorithm 
for adaptive power system stabilisers only to a single machine 
infinite bus model. However, the systems we encounter in 
practice are multimachine systems. The design techniques used 
for single machine model cannot be used for multimachine system 
because of the effect of interconnection between different 

i 

generating units in the latter case. Recently there has been 
much interest in coordinated application and design of PSS in 
multimachine power systems [50~52] . The design for raultima- 
chine PSS is constrained by the fact that only local signals 
are available for stabilization. The control to be used is, 
thus, decentralised control. It has been observed that all 
the machines in a multimachine system need not be equipped 
with PSS [l] . It is enough to provide PSS on a few selected 
units to ensure satisfactory system performance • 

In decentralized control the entire system is 
divided into various subsystems, each having a local conttol 
station. Only local system outputs are used by controllers 
at each control station to control only local inputs. All the 
controller actions are of course, involved in controlling the 



46 


large system. Application of decentralised control has been 
recently done, by Lefebvre [33] using eigenvalue assignment, 
in a 3 machine system with PSS designed on all the machines. 

In the literature, there are algorithms avail- 
able, which consider sequential addition of stabilizers [32,2 
using pole assignment technique. But the sequential additio 
of PSS disturbs the previously assigned eigenvalues. Simulta- 
neous tuning of stabilisers as proposed in Ref. [33] is consi 
dered to be better alternative and used in this work. 

In this chapter we described the system state 
and output equations for a typical, small, multimachine syste 
Coordinated design of PSS is carried out using the technique 
given in [25]. Two different schemes for determining the 
weightages to be applied for adaptive PSS design are presente 
Finally the numerical results are given at a few operating 
points to assess the adaptive PSS design. A schematic block : 
diagram of the adaptive PSS scheme is shown in Fig. 4,1. 

4.2 DFRIGN PROCEDURE 

To start with, we develop the non-adaptive i 
design of PSS based on eigenvalue assignment, as in Anwar’s ; 
thesis [Chapter -6, Ref.l], We then consider the adaptive PSS| 
design. Our design problem for non-adaptive PSS design can ^ 
be stated as: Design PSS for all the machines in the multi-; 

machine system, with decentralised output feedback, to assign 
closed-loop eigenvalues to preassigned locations. 
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Fia. 4,1 Block diagram for Adaptive PSS 
^ in Multimachine Case. 












formulation of control problem 


The state space model of multimachine system 
with n generating units, developed in section 5.2 of Ref. [l. 
can be represented as 

~ ^ ....(4.1) 

where , 

V _ t t it 

—M ^2 ^n ^ 

and 

^ %2 ^mn^ 


represents the state variables of ith generating unit and 
its excitation system. Ay^ is a full matrix and is a block 
diagonal matrix, (see Sec. 5.2, Ref.[l]’'. The state equation 


of the ith generating unit can be written as 


4ii 


1=1 

5^i 


n ^ 2Lr.n- + . U_ • ...(4.2) 

j -mj mi mi 


In eqn.(4.2), and are submatrices of Ay^ correspondin 

to the ith unit and b^^j^ is the ith diagonal block of Simil 

arly the output (eg. rotor speed) of the ith unit, can be 
expressed as 

. 1 

.... (4.3)1 


mi 




Since only decentralised control is considered, the second 
term in eqn. (4.3) can be dropped. The summation term in I 

li 

(4.2) represents the interconnection of the ith generatij 


eqn. 
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unit with other units of the system. The interconnections 
are parameterized with 0< r < 1 as follows 

^ Anii^ 2mi ^ £ [Anij^ 4ij 

A 

Eqn, (4.3/ can be written, in view of decentralised control, 
as 

'^mi ~ ^'^mii^ ^i .,.(4,5) 


with the introduction of parameter r, the degree of inter- 
connection in eqn. (4,4) becomes variable* For r=0, the 
system is completely decoupled and for r=l, it is fully 
connected. 


Initially we design stabilizers for all the machine 
neglecting the interconnection, ( i.e. with r=0 ), with rotor 
velocity feedback, to assign the closed loop eigenvalues using 
the technique of Munro and Hirbod [3l], utilized in chapter 2 
and 3. Then the interconnection is introduced in small steps 
to update the parameters of PSS in such a way that the closed 
loop eigenvalues remain unchanged. Let the state equation and 
control law for PSS for the ith machine be, respectively. 



u, 


mi 


= [Sj] 4 + Ri 

= [Qi] 4 + Kj 


.,,(4.6) 


...(4.7) 


where represents the state variables of PSS«Eqns. (4.4) to 
(4.7) yield the dynamical equation and control law for the 


# 
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augmented system (including PSS dynamics) of the ith unit as 
given below: 



(4.8) 


[C. 1 ] X. 


and 


!il = [K^i ] Yj 


(4.9) 

(4.10) 


where x. = z} ; Y. = [Y . 

—1 -mi —i ■' ' 1 , mi —1 


A. . = 

Anil ° 

* A. .= 

A . . o' 

mij 

11 

_ 0 0 


_ 0 0_ 

Bi 

rbrnl 0~ 


- K, Q- 

■ R. S. 


0 I 


1 1 


''C . . 0 1 

Cii = 1 

I 0 I 1 

For the composite power system, Eqns. (4.8) to (4.10) take 
the form 


X = [Ap + rA^J X + [Bj^] U .....(4.11) 

|. 

Y = [Cj^] X (4.12)j 

1 

U = [Kn(r)]Y (4.23)| 


where Ap = Block diagonal Aini 
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Bq = Block diagonal B^] 

= Block diagonal 

Kjj = Block diagonal [K^j,,K^ 2 ^dn^ 



The objective is to find a control law U , such that the 
composite power system ( eqn. (4.11)), has desired eigenvalues 
The eigenvalues of the composite system are the eigenvalues 
of the closed-loop matrix 

F(r) = Aj 3 -t r/^ + BQKjj(r)CQ 


At r=0, the matrix Kp(0) can be easily computed, since the : 
system is decoupled. Interconnection is then slowly introduc. 
in small steps of Ar, and Kj^(r) is modified such that that i 
the eigenvalues are not much disturbed in the presence of 
interconnections. This is easier than directly computing 


K,,(l), the final answer. i 

i 

The ith eigenvalue equation of the closed loop | 

system can be written, for any value of r as | 


[ A i(r)I-F(r)JV^(r) = 0 



.. (4.14) 



9203 % 
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where Vj^(r) is the right eigenvector associated with the 

eigenvalue A ^(r). The left eigenvector of the system is 
defined by 

l/(r) [Aj(r)l-F(r)] = q .... ( 4 . 15) 

where Wj (r). V^(r) = 1 for j = i 

= 0 for j i 

In the above, * indicates conjugate transpose of the eigenvec— 
^or Wj (r ) . Differentiating eqn. (4.14) with respect to r 
and premultiplying by W*(r) we get 

Wi(r)^ [ Ai(r)l-F(r)] V.(r)+V7*(r) (r )I-F(r ) J V^Cr) =0 

.(4.16) 


It is obvious from eqn. (•,-.15) that the second term of eqn. (4.16 
will be zero. Using the expression for F(r) and rearranging 
eqn (4.16), we get 


IiCbd -ar 




^ - II (4.17) 


In eqn (4.17) 


iiii. 

dr 


may be set equal to zero, so 

dK. 


as to make eigenvalues insensitive to change in r. 

be solved from eqn. (4.17) and utilized in updating as 

follows 


can 


Kp (r+ ..:. r) 


KQ(r) + 


dKQ(r) 

'df 


Ar 


(4.18) i 
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Solution of eqii. (4,17) is obtained by using 
the tensor product operation. Suppose M = is of 

dimension pxq and L is of dimension rxs, then the tensor 
product of M by L (M ® L), is defined to be the prxqs matrix 
[m.jL]. The ’ vec * operation on the matrix M = [Mj^,M 2 » . . . »M^], 
where is the ith column of M, is defined as 


Vec (M) 



M 


q 


It may be shown that the set of equations (4.17) for i=l,2,..,N,. 
where N is the total number of states, is equivalent to 


where 


and 


[0 (r)] Vec [ 


dK2(r) 

dr 


] = . P(r) 


0(r)= 




Q Ml 

[C^V^]^ Q 12'' Bp 


[CpV^]^ Bp 




P(r)= 


J. ^ [A 
d7~ 2 '-AoJVi 


i2 . 12 

dr 


n 

dr 


Wj^[Aq]Vjj 


(4.19) 
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As Kj 3 is a block diagonal matrix, columns of 0(r) corresponding 

dK 


to zero entries of Vec ( -^) are deleted to obtain 0 ( t ) . 
Eqn, ( 4 , 19 ) can now be written as 


0(r) 


Vec ( 
Vec ( 


dK 


dl 


dr 

dK 


d2 


dr 


Vec ( 


dK 


dN 


dr 


) 


= P(r) 


,..(4*20) 


dK. 


From eqn. (4.20), solution for Vec ( ), i=l,...,n, can 


be obtained by computing the pseudo-inverse of matrix 0(r)« 
Solution of eqn. (4.20) in each iteration provides an update 


dK 

of which is used in eqn, (4.18) to update the feedback 

gain matrix Kj^. The conditions for existance of a real 
solution of eqn. (4.20) Vv^re obtained in Ref. [25] to be the 
following 

i) Each subsystem (machine and PSS) is controllable 
and observable. 

ii) The eigenvalues of each subsystem can be assig- 
ned by local dynamic output feedback. 

iii) At most one subsystem has an augmented state 
space which is of odd order. This condition 
guarantees that the feedback gains, i.e. 
matrix K^^, are real. 

iv ) Matrices and have full column and row rank 
respectively. 
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4.4 DESIGN ALGORITHM 


The following algorithm Ref*[lJ was used for 
PSS design in multimachine system. 


Step 1 

Step 2 

Step 3 


Step 4 


Obtain state space equations of each machine as 
in eqns. (4.4) and (4.5). 

Augment state equations of the machines with 
stabiliser dynamics as in eqns. (4.8) and (4.9). 
Design feedback gain matrix each of 

the machines using local velocity feedback, 
to assign closed-loop eigenvalues to desired 
locations, by setting the parameter r=0 in 
eqn. (4.8). 

Construct eqns. (4.11) to (4.13) for the composite 
system from eqns. (4.8) to (4.10) of individual 
machines. 


Step 5 : 


Step 6 ; 


Step 7 ; 


Step 8 ; 


Compute left and right eigenvector of the 
matrix F(r) . 


Calculate 0(r) and p(r) and solve eqn. (4.20) 

for fJS . 
dr 

t 

Update KQ(r) as shown in eqn. (4.18), by 
selecting a sufficiently small value of -i r. 


Increase the value of r to r+ r. 


Step 9 


Go to step 5 if r is less than or equal to 
unity. Print ^^(r) when r > 1 and stop. 
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4 . 5 SYSTEM MODEL 


A multimachine power system problem given in 
Ref. [46], consisting of 3 generators and 9 buses, is considered, 
A single line diagram of the system is shown in Fig, A 4.1. The 
system data used for PSS design are given in Appendix A-4. Each 
machine with its excitation system is represented by a fourth 
order model. Machine 1. is considered as reference for the 
measurements of rotor angle; it, therefore, has a third order 
model. The open loop system model is, thus, of eleventh order, 

4,6 NUMERICAL RESULTS 


We use only velocity feedback to design a second 
order PSS for each machine, with transfer function of the 
form 


© s +0,s+0o 
F(s) = -S i S 


,., (4,21) 


s +Y 1 S+Y 2 

The open loop system matrices are shown in the next page . 
Table 4.1 shows the open loop eigenvalues of the system. 


Following the design technique as refered in 
section 6.4 of Ref [l] for pole assignment neglecting 
interconnections we develop three PSS for the three machines. 
Here we give a description of PSS designed in tabular form. 




•• 0 *; l .293 

0.0 

“ 1 . 8356 

Am . 1 . ; l .™ 

“ 1 . 033.0 

0.0 

0.0 


273 . . 303 

0.0 

“ 100.0 


0 . 0238 

0.0 

“ 0.0063 0.0 


0.9213 

0.0 

12.1999 0.0 


44.236 

0.0 

“ 40.041 0.0 


0.0344 

0.0 

0.0473 0.0 

Am 3 , 3 ™ 

0 . 8548 

0.0 

7.2381 0.0 


65.007 

0 . 0 

75 . 534 0 . 0 

bml ™ 

0.0 0.0 

2500 . 0 

Cm : l . l == 

0.0 1.0 

0.0 



0.01730 

0.0 

0.0 

Am 2 ; l .™ 

■"• 1.7983 

0.0 

0.0 


0.0 

- 1.0 

0.0 


94.2625 

0.0 

0.0 


“ 0.51140 

0.0 

“ 2.65610 “ 1.39120 

Am 22 == 

“ 10.7674 

0.0 

“ 71.0929 0.0 


0.0 

1.0 

0.0 0.0 


215.8138 

0.0 

189.7026 " 66.6667 


0 . 206700 

0.0 

0.00000 0.0 

Am 23 == 

2.6601 

0,0 

22.6567 0.0 


0.0 

0.0 

0.0 0.0 


106.5345 

- 0.0 

90.1385 0.0 

bm 2 = 

0.0 0.0 

0.0 

2500.0 

Cm 22 === 

0.0 1.0 

0.0 

0.0 


0 . 03300 

0.0 

0.0 

Am 3 X -“ 

“ 1.2497 

0.0 

0.0 


0.0 

- 1.0 

0.0 


75 . 6676 

0.0 

0.0 


0.14830 

0.0 

0.997000 0.0 

Am 32 ™ 

5.44670 

0.0 

55.71570 0.0 


0.0 

0.0 

0.0 0.0 


54.6472 

0.0 

“ 125.275 0.0 


“ 0.61980 

0.0 

“ 2.160200 “ 0.9365 

Am 33 " 

“ 26.2893 

0.0 

“ 124.5844 0.0 


0.0 

1.0 

0.0 0.0 


225.4484 

0.0 

371.1237 “ 50.0 

bm 3 ™ 

0.0 0.0 

0.0 

2500.0 

Cm 3 ™ 

0.0 1.0 

0.0 

0.0 



0. EIGENVALUES 

“0 . 427446-1- J8 » 252612 
-0 * 427446 -"J8 . 2526 1 2 
“0 . 880441 -f- J1 1 . 9 1 3736 
”0 ♦ 880441-- J1 1 . 913736 
-~2*068020-fJ0»<) 

--3 * 6S9238-f JO » 0 
~7.976150-f JO.O 
~44»96911-f JO.O 
-61*84220-i-J0.0 
-94*79659-f JO.O 
0*0-fJ0*0 

TABLE 4.1 EIGENVALUES OF THE OPEN LOOP SYSTEM. 


JO. EIGENVALUES 

-0. 26856 1-f JO. 0 
-~1.196420-IJ8. 195561 
-1 . 196420- J8* 195561 
- 1 . 272 1 68-f JO .041 988 
~1.272168--J0. 041988 
~1.593803-iJ12. 150309 
-1 .S93803--J12. 150309 
-2.067331-t-JO.O 
--2. 444427-1 JO. 0 
-3. 035573-1 JO. 0 
-3. 754584-1 JO. 0 
-9. 685784-1 JO. 0 
-16. 377494-1 JO. 0 
-43. 9460 10-1 JO .0 
-48. 069 183-1 JO. 0 
-61 .94 1589-1 JO. 0 
-94. 798347-1 JO. 0 


TABLE 4.5 EIGENVALUES OF THE CLOSED LOOP SYSTEM. 
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Machine 1 

Machine 2 


Machine 3 

Open loop 

0.0 

-0.433541 


-0.983698 

eigenvalues 


+j 8. 230556 


+j 11. 021 860 

at r=0 

-5.393352 

-4.532270 


-3.272960 


-94.735936 

-61.778666 


-45.379431 

Desired 

-2.0, -3.0 

-2.20 +j8.0 


-3.0 + jll.O 

locations 

-4.0 

-3.0 


-3.0 

Remaining 

-1.499999 

-5.840128 


-3.450304 

Closed 

-94.736379 

-13,210880 


-39,758581 

loop 


-48.518170 


-61.888380 

eigenvalues 





Table 4.2 

Eigenvalues 

of the 3-machine 

system at r=0 

Table 4.3 

shows the parameters of the 

three stabilizers 

obtained 

in the first 

stage of design. 

It 

is clear that 

since for 

every machine we specify only 

three locations# 

the pole 

placement is' 

partial . 



Machine 0^ 

Ql 

_ ®2 Ti 


^2 

1. 0.03485 0.121359 

0.0900124 2.885756 

1.456124 

2. 0.299454 1.238136 

0.897746 19.949938 

42.967523 

3. 0.612943 3.417337 

6.138215 51.365486 

151.31207 


Table 4,3 PSS parameters for r = 0. 



60 


The feedback gain matrix for ith machine can be 
obtained in terms of PSS parameters, defined in eqn. (4,21) 
as 



poi 

0 ' 

*=^21 

0 ’ 
^li 



II 

1 — ^ 

0 

0 

1 


' (4.23) 


1 

1 

“Ygi 

""^li 



where 

1 — 

o' - 
^li 

®li 

-®61 

Yii 


and 

II 

•H 

CD 

®2i 

1 

(D 

O 

H* 

^2i 



In the second stage of design, for updating the decen 

tralized feedback gain matrix KQ,a step of 4 r =0,05 is chosen 

to, introduce interconnection. The number of unknown in 

eqn. (4.20) are fifteen The order of closed-loop system 

is seventeen and therefore, eqn. (4.20) represents a set of 

seventeen simultaneous linear algebraic equations in fifteen 

unknowns. Pseudo inverse of 0(r), is determined for obta- 

dK 

ining a solution of in eqn, (4.20). Table 4,4. shows 

the final values of Pss parameters with full interconnection, 
i.e. with r=l. Table 4.5 gives the closed loop eigenvalues 
obtained with PSS designed with interconnection, i.e., with 
PSS parameters given in Table 4.4, 
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Machine 


ei 

®2 

Yl 

^2 

1. 

0.0353845 

0.1178504 

0.0903524 

2.385744 

1.328615 

2. 

0.3969945 

1.3364065 

0.9102773 

19.739929 

43.239750 

3. 

0,6383817 

3.5379695 

6.2815793 

50.264404 

145.92301 


Table 4.4 PSS Parameters for r = 1, 

Table 4,5 gives the closed loop eigenvalues obtained with 
PSS of Table 4.4. 


4.7 DISCUSSIONS 

From Table 4.1 and 4,5 the advantage of pro- 
viding PSS is clear. All the open-loop eigenvalues are pushed 
further into the left half of the S-plane; the oscillatory 
modes, related to rotor oscillation are further damped. The 
algorithm used was computationally efficient and was fairly 
accurate in arithmetic. One point to be noted here is 
that since in designing all the three PSS, we chose second 
order PSS ( unlike in Chapter 2 and 3), for second and third 
machine we had to do partial pole placement instead of 
complete pole placement , It is also observed that the 
remaining solitary unspecified pole of the closed loop system 
got placed well inside the left half plane, so, the performa- 
nce of the closed loop system will not be adversely affected 
by this. However, this entire design was for non-adaptive PSS. 
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In the next section wp shall describe the scheme for adaptive 
PSS, as already outlined in Fig. 4.1. 

4.8 ADAPTIVE PSS FOR MULTIMACHINE SYSTEM 

4.8.1 INTRODUCTION 

As already proposed for the single machine 
case, we design three different PSS for each machine, for 
three distinctly different operating points. In the single 
machine case we had a 2- Dimensional operating space. Even 
in this simple multimachine case, the dimension of the space 
of operating conditions becomes six. So covering the whole 
range of feaseble operating conditions will be a rather 
difficult task and we will need considerably more number of 
PSS designed for each machine to cover the space of operating 
conditions satisfactorily . However, we have not investigated 
into this matter deeply, and our aim in this section is to 
prove that the idea used for design of composite PSS for single 
machine model will also give improved performance, if applied 
for the multimachine model. Two different weighting schemes, 
one based on change in system load and another based on 
machine load as already recommended in chapter 3 are presented 
in the next subsection. 

4.8.2 WEIGHTAGES 

Two different schemes for determining the 
weightages to be given to the different PSS when the operating 
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condition changes are considered. The first scheme is based 

on system loading. As may be observed from table 4,8, among 

the three design points, D~1 is comparatively lightly loaded 

and D “3 is heavily loaded. We have covered a range of 

total real and reactive loading on the system. For any 

system loading, the total complex load power delivered by 

the system is determined. The linear distances of this complex 

power from those at the three design points in the complex power 

plane decides the three weights to be associated with the three 

distances . 

PSS. The weights are inversly propertional to their respective [_ 
Furthermore the sum of the three weightages should equal to 
unity. The numerical results using this proposition are pre- 
sented in Table 4,6, 




w, 

W2 

^^3 

Test Point 1 

System Loads 

Machine 1 

0.4 

0.4 

0.2 

P^ = 120.0, Q^=45.0 


Machine 2 

0.4 

0.4 

0.2 

Pg = 85.0, Qq=25.0 

P^ = 95.0 =30.0 

Machine 3 

0.4 

0.4 

0.2 

Test Point 2 

System Loads 

P^ = 130.0, Q^=55.0 

Machine 1 

0.2 

0.4 

0.4 

Pg = 95.0, Qg = 35.0 

Machine 2 

0.2 

0.4 

0.4 

■P^ = 105.0, Q^= 40.0 

Machine 3 

0,2 

0.4 

0.4 


Table 4.6 V'leights given to the PSS according to scheme 1, 
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The second scheme for determining V\^eights is the 
same as the scheme used in the single machine case. The 
real and reactive powers generated by each machine are found 
out from load flow studies, and for each individual machine, 


the same technique as used 

in section 

3,2, is 

applied, 

Table 

4.7 gives the numerical results using 

this technique at 

the 

same test points. 







*1 

«2 

'^3 

Test Point 1 

Generator 

outputs 

Pgi =56.3728, =^2.185 

Machine 1 

0.4637 

0.4193 

0.1171 

Pg2~163 .0, Qg2 == 6,70 

Pq3= 85.0, Q^3 = -10.9 

Machine 2 

Machine 3 

0.3333 

0.3333 

0.3333 

0.3333 

0.3333 

0.3333 

Test Point 2 

Generator outputs 
Pgi=87.2296,Qgi =61.5441 
P32=163.0, Qq2 =6.70 
%3 =“16).9 

Machine 1 

Machine 2 

Machine 3 

0.2794 

0.3333 

0.3333 

0.4006 

0.3333 

0.3333 

0.3200 

0.3333 

0.3333 


Table 4.7 Weights given to the PSS according to scheme 2, 


4.8.3 ALGORITm FOR ADAPTIVE PSS 

The algorithm used for adaptive PSS can be described 


as follows : 
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Step 1 : Select three different operating points at which 

PSS are to be designed . 

Step 2 : Design PSS for each operating points, following the 

same procedure as in the non adaptive case, 

to 

Step 3 ; Determine weightages£be given to the PSS designed for 

any given operating condition using one of the two 
methods given in previous subsection. 

4.8.4 NUK^ERICAL x^ESULTS 

Following the algorithm described above, numerical 
results have been obtained for two typical operating conditions 
which occur within the range of operating conditions covered by 
the design points considered. Table 4.8 gives the details of 
PSS designed at the three preselected design points. Table 4.9 
gives the details of the composite adaptive PSS for the two test 
points. It also gives a coraparison between use of non-adaptive 
and adaptive PSS. First column in the Table gives the closed 
loop eigenvalues obtained with the non-adaptive PSS designed 
for design point D-l.' Second and third columns of the table 
give the closed loop eigenvalues with the composite adaptive 
PSS, using scheme 1 and scheme 2 respectively, for determining 
the weights. 

4 » 9 • DISCUSSION 

As may be observed from Table 4,9, the improve- 
ment achieved by using the adaptive schemes is not much. It is 
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nevertheless, true that there is some improvement and the 
oscilatory modes related to rotor oscillation are particu- 
larly better damped by the adaptive controller designed. 

It is also true that the number of design points used is 
small. This was because of the complexity involved and 
the time taken for computation. The schemes proposed for 
determining weightages are simple and convenient to use. 
However, they are adhoc and lack theoretical basis. One 
feels that if more number of design points are considered 
so as to cover a wider range of operating condition and a 
better scheme for determining the weightages is used, the 
performance of the closed loop system will show considerable 
improvement. 
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CHAPTER 5 

CONCLUSION 

Traditional nonadaptive PSS designs are satisfactory 
over a small range of operating conditionsaround that for which 
the PSS is designed. There is, therefore, need for developing 
adaptive PSS designs, which hold over a wide range of operating 
conditions. Coventional adaptive designs are complex and 
involve on-line identification of model parameters which undergo 
variation. Practical implementation of such PSS is not simple. 
The need for developing simpler adapation schemes is therefore, 
clear. This thesis has addressed itself to this task. 

The idea used in the design is to have a number of 
PSS, working simultaneously and in parallel. The control 
input to the excitation system is obtained as a weighted sum 
of the outputs of the various PSS. Each PSS is designed for 
a particular operating condition. A wide range of operating 
conditions can be covered by designing a sufficient number of 
PSS. For any given operating condition, the weights for 
obtaining the weighted sum of the outputs. of the various PSS 
are given by a simple rule. 

The above concept of adaptation has been examined in 
detail for a single-machine-infinite-bus system. It has been 
found that a relatively small number of PSS will yield good 
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system performance over a wide range of operating conditions* 
This scheme, is therefore, simple and superior to the traditional 
nonadaptive schemes. It is also better than the traditional 
adaptive schemes from the implementation point of view. The 
additional hardware requirement is not much, in this scheme* 

Design of the proposed adaptive PSS scheme has also 
been examined in a simple multimachine system, involving three 
generators and nine buses. Here, complexity involved in PSS 
design increases significantly. It is also seen that this 
adaptation procedure to choose the relative weights in the PSS 
cannot be chosen, by a simple procedure as in the single 
machine case, with a fair degree of certainty regarding the 
outcome. However, two schemes have been examined, and 
improvement in system performance established. More work needs 
to be done in developing better algorithms for determining the 
weights to be attached to the individual PSS and also for 
determining the number and values of operating conditions 
for which the individual PSS are to be designed. In any case, 
the investigations so far done and reported here appear to be 
promising. 
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APPEI - ^DIX A- 1 

STATE SPACE MODEL FOR SINGLE MACHINE SYSTEM 


Here we consider a three winding synchronous 
machine model. The stator of the machine is represented 
by a dependent current source. A single phase equivalent 
circuit of the power system network is shown in Fig. 2.1 (b). 


The system equations can be written as [34]: 


where 


2H p2n _ j 

Wq me 


T; = (I^i^-I^i^) 


... (A 1.1) 
... (A 1.2) 


Id = Id^/(^'dI'do> •••(A 1.3) 


Iq = ” Iq 

r, = (X'd-Xq)/X'^ 


... (A 1.4) 


From Fig. 2.1 (b): , the armature current plasor i can be 
solved by using the mesh equations 


where 


Let 

where 


A i'^'d ^ A, 

^ “ "(x+xTTT" I “ 


^ “ j ’ (x+x'; ' ) = n + 
i^'d 

aj_+J02 - j (X+X' .) 


Bt — 


(X+X'^)^ 


X' 


(x+x 


D 


...(a l.b) 
... (A 1.6) 
...(A li7) 

... (A 1.8) 
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and 3i2 =0 .... (A 1.9) 

From eqn, ( A 1.5) 

iq-jid=(9i+j32^ sind) .... (A I.IO) 

Equating real and imaginary parts we get 


iq Cosd+y2sind) .... (A l.ll) 

and 

i^ = aj^I^-E(y 2 Ccsd-y^ sind) .... (A 1.12) 


Using eqn. (A 1.4) 

F 

iq= ■ ‘ (y^ Cosd+y 2 sind) .... (A 1.13) 

and 

id = ®i^d“ E(y2^°®^“yiSiiT^^) •••• (A 1.14) 

From eqn. (A 1.16) 

1 

y, = O and y,-, - - , . i 

. I 

In equn. (A 1 J.3) and (A 1.14), i^ and i are expressed in j 

terms of I^ and d only, I 

EXCITATION SYSTEM MODEL 

Many excitation system are used in practice. Ref. [49] | 

gives a detailed survey of all the practical excitation systems. ; 
In this work# a. static excitation system, whose block diagram 
is shown in Fig, 2.1 (c), is considered; 
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In derivation of excitation system dynamics the 
time constant of static exciter may be neglected, since 
it is small. Thus the state equation can be written as 




«*«>( A l.lfj) 


Where 


Vt = [^2 + v 2 ] >/2 


The park's components of the terminal voltage and 


V can be written as [34] 

H 


and 


''d = 


Vq = (Id-id) 


....(a 1.1^] 

... ( A 1.17) 


From Fig 2,1 (c) 


%d ~ "^ref ■" ^1"^^ 


... ( A 1.18) 


LINEARIZED MODEL 


We obtain the linearized model from eqns. (A 1,1), 

(A 1,3) and ( A 1.15). For doing so we eliminate the non-state 
variables Aij and A± which are obtained by linearizing eqns. 

Vj 

(A 1.13), and ( A 1.14). 

^Iq = + KjAS 

A iw = Ko A I . + K. 4 d 


...(A 1.19) 
. . . ( A 1.20) 



ou 


where 

K, = 0 , K^= 

(Vl sind^-y^cosd^), 



l-a3_Ti 




= ajL and = E 

[ - y^ sind^-yj^cosa^J 

... (A 1.21) 

Thus, 

the resulting fourth order system equation can 

be written 

as 





X = AX + bV 


... (A 1.22) 

where 

X = , 

AX, 3" 


and 

V = A u 




A 


*11 


‘31 


*41 


‘12 


0 


‘32 


*42 


0 

1 

0 

0 


*14 

0 

0 

^44 


i b. 


and b= 


I 0 


0 

0 


(A 1.23) 


where 


®11 


) Kg] / 

<^'dT'do) 

®12 

= - Kr(V^'d)/(^'d’''do> 


®14 

= -1/(X 

• T’ ) 
dMo^ 



®31 

“o 

^qo ” 

A 

5<'d 

®32 

■= W (K4lqo-K2ldo-^‘^2i 

doAX'd 

®41 

Kp, 

Tr 

V 

r -22. 
‘ Ao 

XKg ] 


®42 


r Ac 

■Ao 

(E cosd^ 

- XKg) 


^qo^^to 

(-E sind^+XK^)] 

®44 

= - 1/Tf^ 

b 

1 = l/(X’dT'^^) 


(A 1.24) 
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APPENDIX A - 2 

DERIVATION OF SYSTEM MATRICES 

The system data considered are as follows ( .All values are in 
P.U.) 

Generator parameters X^= 1.14, X^ = 0.66, X’^ = 0.24 

T' = 12.0, H = 1.5 , m = 314.0 
do ' o 

Kj^ = 50.0, T^= 0.01 sec. 

Network Parameters R=0.0,X=0.7 

Operating Data E = 1.0 ^ 0° 

P = Received real power = 0.8 
Q = Received reactive power =0.6 

The elements of the system matrices can be calculated in the 
following way. 

+ Zi 

= 0.58+j0.56 = 0.8062 Z 43.995 

Therefore 6^= 43.995° 

Now , Pg + jQg = 0Q . i^ 

Therefore p = 0.8 and Q = 0.1 

y y 


= 0.8 + jO.l 
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Furthermore , 


Q. 


0^ = tan ^ ) = 7.125 


o 


’ao 


Pg/Vt COS 0Q = 1.0 = i 


T i^^X^cos0^--i^„r_ sin0^ 

B = tan-1 ., a° _ q .._.. ':o ap _a o j. 


= 36.407^ 


= Pgo+®o=®0-^02° 


''do = ''t ^1" Pao = °-^285 
r 

qo 


go 

v^.. =V^ cos = 0.6488 


Ido = lao =l"(Pgo-«'*o^ = 

Iqo = lao (Pgo+e>o) = °-''250 


V 

I __go 

do 

^ d 


+ = 3.3921 


Iqo =1 Iqo =- 1- 26825 


Now using the formulae given in Appendix A-1 we get 


11 

=-0.16312 

®12 ” 0.3278 

®14 

31 

=-26.34 

3^2 = “ 47.5885 

®41 

■42 

=-773.2739 

^44 ~ ”* tOO.O 

^1 


= -0.34722 
= 719.0961 


0.34722 



APPENDIX A - 3 


POLE ASSIGNMENT WITH OUTPUT FEEDBACK f3l' 


The technique given here is based on the work of 
Munro and Hirbod [31 j. The method presented in [31 ] provides 
a technique for the design of full-rank compensators for multi- 
varible linear systems. Here the technique is restricted to 
single input linear systems 


and 


X = [aJ X + b 


Y = [Cj X 


u 


... (A 3.1) 


... (A 3.2) 


where, X£ R^, u € and Y£ R^ are respectively vectors of 
state, input and output variables. The r X 1 open-loop trans- 
fer function matrix G(s) is 


G(s) = C(sl-A)"^ b= 

^o'‘ 




Nj_(s) 


N^(s) 


... (A 3,3) 


where 


[J-qCs) =s + ... +a^ 


and 


^i^®^ - ^il®"”^ +Pi2®^""^ + ••• + Pin 


....(A, 3.4) 


.... (a 3,5) 


lJ,^(s) is the characteristic polynomial of G(s), Consider the 
output feedback law 
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where R'^ represents the reference input . 

The problem can be defined as the determination of 
1 X r dynamic feedback compensator, F(s), such that the 
closed-loop system 

H(s) = [I + G(s) F(s)r^ G(s) .... (A 3.7) 

has a desired set of eigenvalues. The kth order compensator 
F(s) has the form 

F(s) = [Mj^(s),M 2 (s), ..., M^(s)J ...(A 3.8) 

where 

Pj,(s) = s +Y2S + Yg® ^ ••• Y|j ...(a 3 , 9 ) 

and 

M^(s) = ®ik (A 3.10) 


u (s) is the characteristic polynomial of compensator transfer 
function F(s) . 


The resulting closed -loop system characteristic polynomial 
M'^(s) , defined as : 




n+k 


djS 


n+k-1 


,+ d 


n+k 


(A 3.11) 


can be written as [31 ] 

= P'o(s). 


r 

tx^(s) + S N. (s).M. (s) 
i=l ^ ^ 


(A 3.12) 
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COMPLETE POLE ASSIGNMENT 

The problem of pole assignment can be defined as: 

given |j,^(s), N^Cs), N^(s) and a desired set of closed- 

loop poles, find p (s) , M-,(s), M (s) of lowest degree 

which satisfy eqn. (A 3.12). Equating coefficients of like 


powers 

of 

s in eqn. ( 

A 3.12), 

we get 






0 ^ 

1 — 1 

1 I 

= 



* • • » 

. (A 3.13) 

where 









1 

0 ... 0 

^11 •• 

. 0 

^rl 

...0 

i 

1 


“l 

1 ... 0 

Pl2 •• 

. 0 

^r2 

. . ,0 

1 


“l 

“l 


• 





* 

“2 ^ 


Pll 


Prl 

1 

II 

# 

* CX 

^In 

^12 

•^rn 

^r2 



°^n 

• ^*2 

0 


0 




0 

“n 

0 


0 




0 

• 

0 







• 

0 

0 oc 

n 

0 

^in 

0 

Prni 



— 





— J 



k - 

columns 

k+1 columns 


k+1 columns 


Pr 

^ ^1 ^2 ** 

<D 

O 

. .. 0 

Ikl • 

• • i 

' ro 

0 1^ 
••• ^rk-^ 

and 










II 

{ 1 

1 

M 

Csl 

1 

04 


i(d„- 

“n^l ‘^n+l*“^n+k^ 
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Equation (A 3.13) is a set of (n+k) equations in [(k+l)r+kj 
unknown parameters of F(s). The difference vector contains 
the coefficients of the polynomial 

= p^(s) - Pq(s). s^ .... ( A3. 14) 

A necessary and sufficient condition for the existence of 
solution of eqn. ( A 3.13) is 

Rank [Xj^] = Rank [Xj^, .... (A 3,15) 

and a solution of eqn, ( A3*13) is 

4 = C Xg^] dj^ .... (A 3.16) 

gi 

where ^ is the generalized inverse of Xj^ such that 

[xj [X9^] [X^] = [x^] 

PARTIAL POLE ASSIGMV4ENT 

Consider the case when q poles are to be assigned, 
t = (n+k-q) poles are allowed to assume arbitrary locations. 
For this case the closed-loop system characteristic 
polynomial p^(s) has the form 

[-1^ ( s ) = ( s*^+dj^ s^ +. * » + d^ ) ( s +e j^s +...'^ 

.... (A 3.17) 

where •**» *^q specified and e 2 ,e 2 » are to 

be determined. Here the difference vector dj^ can be obtained 
by equating the coefficients of like powers of s in (A 3.17) 
and (a 3,12) as 
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~k ^ -1 (A 3,18) 

where the vectors and ( i=l, t ) contain res- 

pectively the coefficients of the polynomials 


and 


Eqn, ( 


where 


and 


’ Cs) = |i^(s).s^ -p^(s).s^ (A 3.19) 

l-i'i'(s) = p (s).s^'"^, (i = l, t) .... (A3. 20) 

X M 

A 3.13 ) in this case will take the form 

[X,^] Pk = d + [Dj e (A 3.21) 

[d] = [ ^2 ^-^2^ 


Eqn. (A 3.21) can be rearranged as 

= 

where 

[X\] = [Xj^,“D] 

£i; = [pj . ( A 3.23) 

Eqn. (A 3.22) is a set of (n+k) simultaneous algebraic equations 
in [r(k+l) + k+tj unknowns. Vector P'j^ contains parameters of 
F(s) and the coefficients of polynomial of unassigned closed- 
loop poles. The necessary and sufficient condition for the 
existence of solution of eqn. ( A 3.22) is given in eqn (A 3.15) 
with Xj^jPj^ and replaced by and d*j^ respectively. 
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APPENDIX A-4 

3 A4ACHINE 9~BUS POWER SYSTEM DATA 

The single-line diagram of a three machine nine 
bus typical power system is shown in Fig. A4.1, Table A4.1 
gives the line data. Generator and AVR data are presented 
in Table A4.2. Data given in these two Tables are in per 
unit on a base of 100 MVA. 


Line 

No. 

From 

bus 

To 

bus 

Line Line 

Resistance Reactance 
(R) (X) 

Half line 

charging 

Admittance 

Turns 

Ratio 

1 

1 

4 

0.0 

0.0576 

- 

1.0 

2 

2 

7 

0.0 

0.0625 

- 

1.0 

3 

3 

9 

0.0 

0.0586 


1.0 

4 

4 

5 

0.0100 

0.850 

0.0880 

— 

5 

5 

7 

0.0320 

0.1610 

0.1530 

- 

6 

6 

9 

0.0390 

0.1700 

0.1790 

- 

7 

7 

8 

0.0085 

0.0720 

0.0745 


8 

8 

9 

0.0119 

0.1008 

0.1045 


9 

4 

6 

0.0170 

0.0920 

0.0790 

- 


Table A 4.1 Line and transformer data for 3-machine 
system model. 
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Generator 

1 

2 

3 

Rated MVA 

245.5 

192.0 

128.0 

KV 

16.5 

18.0 

13.8 

Power factor 

1.0 

0.85 

0.85 

Type 

Hydro 

Steam 

Steam 

Speed (rpm) 

1800 

3600 

3600 

X^(P.U. ) 

0.1460 

0.8958 

1.3125 

X’^CP.U.) 

0.0608 

0.1198 

0.1813 

X„(P.U. ) 

0.0969 

0.8645 

1.2578 

q 




X^(P.U.) 

0.0969 

0.1969 

0.2500 

leakage (P.U.) 

0.0 

0.0 

0.0 

T'do (sec) 

8.96 

6.00 

■ 5.89 

T'qo 

0.0 

0.535 

0.600 

H(MWS/100 MVA) Sec 

23.64 

6.40 

3.01 

Vj^ Gain 

50.0 

50.0 

50.0 

Tp^(Sec) 

0.010 

0.015 

0.020 


Table A 4.2 . Generator and voltage regulator data for 
3- machine system model. 
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